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We consider P systems with a linear membrane structure working on objects over a unary alphabet 
using sets of rules resembling homomorphisms. Such a restricted variant of P systems allows for 
a unique minimal representation of the generated unary language and in that way for an effective 
solution of the equivalence problem. Moreover, we examine the descriptional complexity of unary P 
systems with respect to the number of membranes. 

1 Introduction 

One decade ago, inspired by biological systems to be found in nature, Gheorghe Paun introduced mem- 
brane systems (e. g., see [4]) as a class of distributed parallel computing models, originally working on 
multisets of objects. The main feature of membrane systems (soon called P systems in the literature) 
are the hierarchically structured membranes embedded in the outermost skin membrane. Every mem- 
brane encloses a region possibly containing other membranes as well as multisets of specific objects that 
evolve according to evolution rules associated with the regions. A computation is a sequence of con- 
secutive configurations of the system evolving by applying the evolution rules in parallel to the objects 
in every region in a maximal manner. The vector of the multiplicities of objects present in the specified 
output membrane in the final configuration of a halting computation is considered to be the result of such 
a computation; in that way, a membrane system computes a set of vectors of natural numbers. 

Many variants of P systems are able to characterize the family of recursively enumerable sets of 
vectors of natural numbers (the Parikh sets associated with recursively enumerable languages). The first 
detailed definitions and results can be found in [4], and an overview of many variants is given in [6]. The 
first monograph on membrane systems @ shows the great variety and potentials of P systems; for the 
actual state of the art in membrane systems we refer the interested reader to the P page O. 

In Q it is shown that special variants of P systems with the underlying membrane structure being 
only a linear tree can be computationally complete when using adequate evolution rules. In 0, several 
quite simple variants of such P systems were investigated that allowed for establishing an infinite hier- 
archy with respect to the number of membranes, which usually cannot be the case for computationally 
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complete variants of P systems. These special variants of P systems were inspired by the tree systems of 
morphisms (introduced in O) and the compound Lindenmayer systems (investigated in (H). In this pa- 
per, we consider those variants as introduced in (3j , yet especially over a unary alphabet, and investigate 
the descriptional complexity of the unary languages generated by such systems. 

The rest of the paper is organized as follows: In the second section we give some preliminary defi- 
nitions of notions needed in this paper and then define the special variant of self-reproducing P systems 
to be investigated in the succeeding sections; moreover, we also construct a suitable example and prove 
some first results. In the third section, we show how to construct a unique minimal representation of 
the unary language generated by a unary P system, which also allows us to show the decidability of the 
equivalence problem for unary P systems in the succeeding section. The descriptional complexity of 
unary P systems with respect to the number of membranes is investigated in the fifth section. 

2 A specific variant of P systems 

In this section, after having specified some notions, we introduce the main variant of P systems investi- 
gated in this paper, i. e., P systems with a linear membrane structure and rules mapping each object to 
powers of itself. Moreover, we give an illustrative example and prove two useful lemmas valid for this 
specific variant of P systems. 

2.1 Preliminaries 

For the basic notions of formal language theory we refer the reader to monographs such as JH. We just 
mention some specific notions used throughout the paper: The empty string is denoted by A, and for the 
length of a string w we write \w\. The set of mappings from a set M to a set N is denoted by N M . 
We denote the set of positive integers by N + and the set of non-negative integers by N; N ra denotes the 
set of vectors (x\, . . . ,x n ) with xi £ N, 1 < i < n, and for (x\, . . . ,x n ) G N ra , (x\ , . . . , x n ) T denotes the 
corresponding transposed vector. Moreover, we write P for the set of prime numbers, and throughout the 
paper, pi denotes the i-th prime number, i. e., p\ = 2, pi = 3, pj = 5, . . . 

2.2 Self-reproducing P systems 

Now let us consider the following very special variant of P systems, where 

• the underlying membrane structure is a linear tree [o[i - • • [n ]n • • - ]i]o; 

• the only axiom is in the innermost membrane; 

• except for the skin membrane (labeled by 0), each region i > 1 contains the membrane dissolving 
operation S as well as for each object a exactly one (noncooperative) rule of the special form 

a — > a m (* ,a ) for some m(i,a) > 0; 

• there are no rules in the skin membrane; 

• the result of a computation is the multiset finally appearing in the skin membrane (which consti- 
tutes a halting computation). 

There are no additional features used when applying the (noncooperative) rules of the special form 
a — y a m ( J,a ) in a maximally parallel manner; for example, there are no priority relations among the rules. 
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The movement of objects from the innermost membrane to the skin membrane only happens by 
dissolving one membrane after the other one using the membrane dissolving operation 5 which may be 
chosen at any time for a membrane containing at least one object. 

These special P systems as described above were introduced in Q and had been inspired by the 
tree systems of morphisms (introduced in [2]) and the compound Lindenmayer systems (investigated 
in HI). In the case of linear trees, these tree systems with morphisms intuitively correspond with the 
self-reproducing deterministic linear P systems as defined above. The main difference is the (additional) 
skin membrane, which only serves for terminating the computation and collecting its result. 

Due to this correspondence with tree systems with morphisms, we can easily describe a P system of 
the special form as defined above in the following way: 

n = (V,hi,...,h n ,w) 

where V is an alphabet, n is the height of the linear tree describing the membrane structure, w is the 
axiom over V, and the hi are the homomorphisms defined by the rules of the form a — > a m( - l ' a ^ in region 
i. As each object is mapped to powers of itself only, we call these P systems self-reproducing. 

As it is easy to see, the language generated by a self-reproducing P system 

n = (y,hi,...,h n ,w) 

can be written as 
in other words, 

L,(n) = {/ l ™ l (fer-(CW)-)|m t >0,l<Kfi}. 

We may also demand that every homomorphism has to be applied at least once before the correspond- 
ing membrane may be dissolved; in that case, the corresponding language generated by 
n = (V,h\, . . . ,h n ,w) can be written as 

ht(h+--{h+(w))--), 

in other words, 

L+(n) = {ti? { (h™ 1 ■ ■ • (/C(«0) ■ ••) I m i > M < » < n )- 

2.3 An illustrative example and first results 
Example 1. For each n > 1, consider the language 

L (n) = {a\ l ...a p n n \ rm > 0, 1 < i < n} 

where pi denotes the z-th prime number. 
L (n) is generated by the P system 

n(n) = (V,hi,...,h n ,w) 
with V = {ai, . . . , a n } , w = a\ . . . a n , hi(ai) = a Pi and hi(ak) = a-k for k / i, 1 < k < n, 1 < % < n. 

m l m n 

Obviously, a\ l ...a£ n " is obtained as /i™ 1 (••• {h™ n (ai...a n )) ■■■ ); moreover, only elements from 
L (n) can be generated by n (n), hence, L (n) = L* (n (n)). Finally, we observe that 

L+(II(n)) = {a\' ...al n \mi> 1,1 <i<n}. o 
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As it was shown in [3], L(n) cannot be generated by a self -reproducing P system with less than n 
homomorphisms (or equivalently, less than n + 1 membranes), i. e., these variants of P systems establish 
an infinite hierarchy with respect to the number of membranes. 

We now establish two simple but nevertheless useful lemmas. The first observation is that the order- 
ing of the homomorphisms is irrelevant (which directly follows from the commutativity of multiplica- 
tion): 

Lemma 2. Let IT = (V, h\ , . . . , h n , w) be a self-reproducing P system, u be a permutation of 1 , . . . ,n, 
and IT = (V, h u n\,. . . , h u ^,w). Then -L*(IT) = L*(IT) as well as L + (IY) = L + (TT). 

Proof: Let a € V and let |u;| a denote the number of objects a in the axiom w of IT; then obviously 

h u(l) ] (C(2) 2) • ■ • (Cfnf' ( aK 0) • • = h T (h? 2 ■ ■ ■ (C («'"'•)) • • • ), 

as, due to the commutativity of multiplication, (a m ) k = a mk = a km = (a k ) m for all k,rn > 1. □ 

Using Lemma [U next we show that we may move the first application of every homomorphism into 
the axiom, i. e., for every self-reproducing P system IT there is an equivalent self-reproducing P system IT 
such that = L*(IT). 

Lemma 3. Let II be a self-reproducing P system. Then we can effectively construct a self-reproducing P 
system IT such that L + (IT) = L*(IT). 

Proof: For IT = (V,h\, . . . ,h n ,w) we obtain 

L + (IT) = {h™ l (h™ 2 ---(h^(w))---)\m i >l,l<i<n} 

= {h^- 1 (h^- 1 ■ ■ ■ (C"- 1 (h (• • • (K(w)) •••)))■••) I m > 1, 1 < * < n} 
= {hT' 1 (hf ■ ■ ■ (C" («/)) • • • ) I < > 0, 1 < i < n} 

= L,(n'), 

where IT = (V,/ti, .. .,h n ,w') with w' = h\(- ■ ■ (h n (w)) ■■■). □ 
The following example shows that the converse of Lemma [3] does not hold: 

Example 4. Consider self -reproducing P systems over a one letter alphabet, i. e., unary P systems of 
the form IT = ({a} , h\ , . . . , h n , w) as considered in more detail in the succeeding section; without loss of 
generality, let us assume that none of the homomorphisms hi with 1 < i < n equals the identity hi (a) = a, 
because it would have no effect on the results. Clearly, the axiom w belongs to the language L*(n). 
Moreover, there cannot be any other string in L*(IT) that is shorter than the axiom. So, choosing w = a Pi 
and a sole homomorphism h\ such that h\(a) = a Pj , where pi and pj are different prime numbers, we 
obtain a system IT = ({a} , h\ , a Pl ) whose language L* (IT) is {a PiP J j m > 0}. On the other hand, it is 
easy to see that for any other unary system IT the string a Pi cannot belong to L + (IT) unless the axiom is 
a and in IT there must be a homomorphism h! with h! (a) = a Pi . But in this case, a Pi belongs to L + (IT), 

2 

too. This is a contradiction, since pi and pj are different prime numbers and, thus, a p * cannot belong 
toL*(IT). o 
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We now consider unary P systems, i. e., systems of the form 



IT= ({a},hi,...,h n ,w) 



where n > 1, w G {a}* is the axiom, and hi : {a} — > {a}*, 1 < % < n, are homomorphisms. Depending 
on whether or not each homomorphism has to be applied at least once, as above we distinguish two 
languages generated by II: 



Now we turn to a language representation which is advantageous for our purposes. To this end, we 
consider the decomposition of positive integers into prime factors. We denote by F the set of mappings 
from P to N with finite support, i. e., only finitely many primes are mapped to non-zero values. Then the 
decompositions are represented by the mapping pf : N + — > F. For example, let n = q\ l ■ ■ ■ q e k k , where 
k > 1, qi prime, qi / qj for i ^ j, ej > 1, 1 < i, j < k. Then {pf{n))(qi) = e^, for 1 < i < k, and 
(pf(n))(q) = 0, for q £ {q\, . . . For convenience, we write pf n for the mapping pf(n). Clearly, 
it holds pfmn = pfm+pfn- Next, let gpf{m) be the index of the greatest prime factor with non-zero 
exponent in the decomposition of m, if m > 1, and set gpf(l) = 1. 

Since pf is injective, every unary language L is uniquely represented by the set 



(Recall that pi always denotes the z-th prime number.) Alternatively, V may be a set of vectors of natural 
numbers, e. g., (mi, . . . ,m n ) T G V implies a m E P~ l {V), where m = p™ 1 ■ ■■p™ n - From this point of 
view, a vector of natural numbers may be seen as a mapping from F and vice versa. In the following, we 
will use both notions in a synonymous way. 

Let IT = ({a},w,h\, . . . ,h n ) be a system, where hi{a) = a c \ for a > 1, 1 < i < n. We set 
k = max{gpf(\w\),gpf(ci), . . . gpf(c n )} to be the largest prime factor with non-zero exponent appear- 
ing in | w | or one of the exponents q, and define the matrix: 



£*(n) 
L + (n) 



{^(^••■(C"H)---) \m,i>0,l <i<n}, 
{^'(^■■•(CW)-) I mi > 1,1 <i<n}. 



P(L) = {pf m \a m GL}. 



Conversely, since pf is surjective, every subset V Q F uniquely represents a unary language 



p-\V) = {a m | there exists / G V such that m = p{ [n) p{ [n) p{ {n) ■■■}. 



( P/ci(Pl) PfdiPl) Pfcn(Pl)\ 
PfrAPl) Pfc 2 {P2) ■■■ PfcSPl) 



\Pfa(Pk) Pfc 2 (Pk) ■■■ PfcSPk) j 



Lemma 5. Let IT = ({a}, h\, . . . ,h n ,w) be a unary P system. Then 

L + {U) = p-\M n {n n + f + b) and L*(IT) = p-\M n {W l f + b) 
where b = (pf\ w \{pi),pf\ w \(p 2 ), ••• ,pf\ w \{Pk)) 7 
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Proof: Let k = max{gpf(\w\),gpf(ci),. . -,gpf(c n )} and foj(a) = a c \ for Cj > 1, 1 < i < n, be defined 
as above. There exist {m\ m n ) € and {e\ , . . . , e^) T € Mn(N™ ) T + b such that the following is 
true: 

a m £L+{U) 



a 



m -h^(h^--(h^(w))--) 



m = |to|c„ C n _j ■ • -Cj 



Pfm(Pi) = Pf\w\ (.Pi) + m nPfc n (Pi) H h mip/ Cl (pj), 1 < i < fe 

pfmiPi) = (M u (m u m 2 ,...,m n ) T + b)\i], I <i<k 

Pfm{Pi)=ei, l<i<k, pf m {q)=0, q <£ {pi, . . . ,p k } 

ei fit 

<^=^ m=p l ■ ■ -p k 

a m eP- l (M n (Nl) T + b). 

The assertion for L*(TI) follows in an analogous way. □ 
Due to the pumping lemma for regular languages, we immediately observe the following result: 

Corollary 6. The only context-free languages that can be generated by a unary P system are of the 
form {a m }for some m > 1. 



4 Equivalence and minimality of unary P systems 

We now may take advantage of the chosen language representation in order to show the decidability of 
the equivalence problem for unary P systems. Moreover, we obtain a unique minimal representation. 
Consider the partial ordering 

f<9^ f{p) < 9{p) for ah> € P 

defined on the set of mappings N p . As usual, we call an element / minimal if g < f implies the iden- 
tity g = f. Let IT = ({a} , hi, . . . , h n ,w) be a unary P system and V = M n (N n ) T + b be such that 

(IT) = P~ l (V). In any other system IT = ({a} ,h[,..., h' n ,w') which generates the same language, 
i. e., V = M n / (N n ') T + b' such that P" 1 (V^ 1 ) = (V), the same prime factors as in n must occur. 
Hence, Mn and M n ' must have the same number of rows. Moreover, a decomposition b = b\ + 62 implies 
b > b\. Therefore, b is the unique minimal element in V and, thus, the unique minimal element in V', 
which implies b = b' . For this reason, in order to deal with equivalence and minimality, for what follows 
we may assume 6 = and consequently may minimize and compare sets of the form Mn(N ra ) T . 

An element x G Mn(N ra ) T \ {0^} - where 0& denotes the zero- vector from (N ra ) T - is called irre- 
ducible, if there is no decomposition of the form x = X1+X2, where x\,X2 £ Mn (N n ) T \ {0& }. The next 
theorem leads to a unique minimal representation. 

Theorem 7. Let Ube a unary P system with n homomorphisms. If M' is the matrix whose columns are 
the n' irreducible elements o/Mn(N n ) T \ {0^}, then 

(i) all columns of M' are columns of Mn and 

(ii) M / (N n ') T = M n (N n ) T . 
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Proof: Since an irreducible element of Mn(N n ) T \ {0} cannot be decomposed into summands, it has a 
representation of the form Mn(0, . . . ,0, 1,0, . . . ,0) T . Therefore, it is a column of Mn, and assertion (i) 
follows. 

By (i), the inclusion M'(N n ) T C Mn(N n ) T follows immediately. Now assume that in contrast to 
assertion (ii), the inclusion is a proper one. Then there is a minimal element x G M n (N n ) T \M / (N ra ') T . 
Since 0& G Mn(N n ) T flM'(N n ) T , we have x / 0&. Moreover, x is not irreducible, since otherwise 
it would be a column of M' and, thus, would belong to M'(N n ) T . So, there is a decomposition 
x = x\ + X2, where x\ , X2 G Mn(N n ) T \ {0}. This implies x\ < x as well as X2 < x. Since x is minimal 
in M n (n n ) T \M'(N n ') T , both elements x x and x 2 belong to M'{W') T . Obviously, M'(W') T is closed 
under addition. Therefore, x belongs to M'(N n ) T , too. From this contradiction we infer assertion (ii). 

□ 

The preceding theorem shows that M'(N n ') T + b is the unique minimal representation of the lan- 
guage generated by IT, i. e., L*(IT) = P~ l (M'(N n ) T + b). It is easy to determine the irreducible columns 
out of the n columns of Mn. Hence, matrix M' can effectively be constructed from Mn. Conversely, 
given some M' with n! columns and b, the unary P system IT' whose language is P _1 (M'(N n ) T + b) 
can effectively be constructed, too. 

Corollary 8. There exists an effective algorithm which minimizes the number of homomorphisms (mem- 
branes) of a unary P system. 

Corollary 9. The equivalence of unary P systems is decidable. 

5 Descriptional complexity of unary P systems 

This section is devoted to descriptional complexity issues of unary P systems with a different number 
of homomorphisms (membranes). The key question is how succinct a language, given by some unary P 
system with n homomorphisms, can be represented by some unary P system with at least n homomor- 
phisms. 

In order to talk about the economy of descriptions we first have to define what is meant by the size 
of a system. In general, we are interested to measure the length of the string that defines a system. In 
particular, we use more convenient size measures, such that there is a recursive upper bound for the 
length of the defining string depending on the chosen size measure. For example, the size of a finite 
automaton equals the product of the number of its states and the number of its input symbols. 

The size |IT| of a unary P system IT = ({a} , hi, . . . , h n , w) is defined to be 

|n| = M+ £ \hi(a)\. 

l<i<n 

We denote the family of unary P systems by T. Clearly, the considered measure implies a total, 
recursive function mapping a unary P system to its size, such that T is recursively enumerable in order 
of increasing size, and does not contain infinitely many members of the same size. 

Let T\ and Ti be two subfamilies of T . A function / : N + — > N + , with fin) > n, is said to be an 
upper bound for the increase in size when changing from a minimal description in T\ to an equivalent 
minimal description in Ti, if 

min{|IT| | n G Ti generates L} < /(min{|IT| | n G T\ generates L}) 

for all languages L generated by some system in T\ as well as by some system in T%. 
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The following theorem is an immediate consequence of Theorem |7J Since any matrix of some rep- 
resentation contains at least the columns of irreducible elements, any equivalent system has at least the 
homomorphisms associated with these columns, respectively. 

Theorem 10. Let Ylbe a unary P system with n homomorphisms where 

L4U)=p-\M n (n n ) T + b). 

If the columns of Mu are different irreducible elements of Mn(N n ) T \ {0&}, then every unary system IT 
with L*(n) = L*(rT) is at most of size \Yl\. 

So, in case of L*(IT), the identity is an upper bound for the increase of size when possibly changing 
from a description with n homomorphisms to an equivalent description with n! < n homomorphisms. 
The situation is different in case of L + (IT) : 

Theorem 11. Let Ylbe a unary P system with n homomorphisms. If there is a unary P system IT' with 
n — 1 homomorphisms such that L + (TT) = L + (lT), then |TT'| £ 0(jTT| 2 ). 

Proof: If a language L = L + (Yl) is generated by some unary P system 

IT= ({a},hi,...,h n ,w) 

whose number of homomorphisms can be reduced, then there is at least one column of M n which is not 
an irreducible element of Afn(N n ) T \ {0^}. Let h n be the associated homomorphism, then 

IT = ({a} ,hi, . . . , /i n _i , h n (w)) 

generates L, too, i. e., L = L + (IT). Therefore, we can approximate the size of IT' as follows: 

Let m = max{|/ij(a)| | 1 < i < n}, then 





|TT|- 


\K(a)\ 


— \w\ + \h n (w 


< 


|IT| - 


\K(a)\ 


— \w\+m \w\ 


< 


|IT| - 


1 + (m 


-i)H 


< 


|IT| - 


i + (|n 


-i)*|n| 




|IT| 2 - 


-ieO( 


|n| 2 ). 



The next lemma shows that the upper bound of Theorem QT] can be reached. 

Lemma 12. Let IT = ({a} , hi , h2,w) where h\ (a) = /12(a) = a m and w = a m , for some m > 2, then 
any unary P system IT with L + (TT) = L + (lT) has size n(|TT| 2 ). 

Proof: Obviously, h\ (a) = a m is associated with a necessary column in Mu- So, the sole equivalent 
unary P system with one homomoiphism is 

W = ({a},h u h 2 (w)). 

For the sizes we obtain |IT| = 3m and |TT| = m + m 2 € n([TT| 2 ). □ 
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If we reduce the number of homomorphisms by more than one, i. e., we iterate the construction x 
times, then we do not obtain sizes of the form Odnl 2 *). The reason is that we have to multiply the sizes 
by at most a factor | O | in each step. 

Corollary 13. Let Yl be a unary P system with n homomorphisms. If there is a unary P system II' with 
n — x homomorphisms such that L + (II) = L + (JV), then |IT| G 0(|II| :c+1 ). 

By an immediate generalization of the proof of Lemma[l2]we obtain a matching lower bound (in the 
order of magnitude) in the worst case. 

Lemma 14. LetH= ({a} , hi, . ■ . ,h n ,w) where h\ (a) = ■ • • = h n {a) = a m andw = a m , for some m>2, 
then any unary P system IT with n — x homomorphisms, 1 < x < n — 1, and L + (JY) = L + (Yl') has 
size a(\U\ x+1 ). 

Finally, if we fix the number m and consider arbitrarily large n, then the trade-off in the size may 
become exponential with respect to the base m. 

Example 15. Let IT = ({a} , h\ , . . . , h n , w) where h\(a) = ■ ■ ■ = h n (a) = a m and w = a m , for some 
m>2, then any unary P system IT with one homomorphism and L + (H) = L + (IT) has size ndnl"). o 

The representation and the descriptional complexity of languages generated by less restricted variants 
of P systems remains as a challenging task for future research; for example, we may ask to which extent 
the results proved above for unary P systems can also be formulated for self-reproducing P systems. 
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